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1 INTRODUCTION 

> 

. ^ ^ The central accretors in a large number of low mass X-ray 
■ binaries (LMXB) are believed to be neutron stars, rotating 
$_( ' rapidly due to accretion-induced angular momentum trans- 
fer. LMXBs are thought to be the progenitors of milli-second 
(ms) radio pulsars (Bhattacharya & van den Heuvel 1991) 
like PSR 1937-^21 with P ~ 1.56 ms (Backer et al 1982). 
The recent discovery of ms {P ~ 2.49 ms) X-ray pulsa- 
tions in XTE J1808-369 (Wijnands & van der Klis 1998) 
has strengthened this hypothesis. Kilohertz Quasi-Periodic 
Oscillations (kHz QPOs) seen in a number of LMXBs is an- 
other indicator of the rapid rotation of the accreting neutron 
star. For example, in beat-frequency models of kHz QPO, 
the difference (~ 300 — 500 Hz) in the frequencies of two 
simultaneously observed kHz QPO peaks is interpreted as 
the rotational frequency of the neutron star (e.g. van der 
Klis 2000). 

Fast rotation makes the neutron star considerably 
oblate and also modifies, both the interior and exterior 
space-time geometry. These, in turn, modify the size and 
the temperature profile of the accretion disc and hence the 
corresponding spectrum (Bhattacharyya et al. 2000;2001). 
General relativity also plays an important role in shaping 
the spectrum. As neutron stars are very compact objects, 
the effect of general relativity is very significant, particu- 
larly near the surface of the star. 

For luminous LMXBs, the observed X-ray spectrum can 
be well fitted by the sum of a multi-colour blackbody spec- 
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trum (presumably from the accretion disc) and a single tem- 
perature blackbody spectrum (presumably from the bound- 
ary layer) (see Mitsuda et al. 1984). The multi-colour spec- 
trum can be calculated if the temperature profile in the ac- 
cretion disc is known. Recently Bhattacharyya et al. (2000) 
have calculated the disc temperature profile as a function of 
the spin rate (fi) and the mass (M) of the neutron star for 
different equations of state (EOS), including the effects of 
rotation and general relativity. In this paper, we calculate 
the observed spectrum from the accretion disc, using the 
temperature profiles. While doing so, we take into account 
the effect of gravitational redshift, Doppler broadening and 
the light-bending effect in the gravitational field of the ac- 
creting star. 

For non-rotating neutron stars, an attempt has been 
made by Ebisawa et al. (1991) to compute the disc spec- 
trum including the effects mentioned above. Sun & Malkan 
(1989) included relativistic effects of disc inclination, includ- 
ing Doppler boosting, gravitational focussing, and gravita- 
tional redshift, on the observed disc spectra for both Kerr 
and Schwarzschild black holes (in the context of supermas- 
sive black holes). They find (as we also do) that the higher 
energy part of the spectrum gets significantly modified due 
to the general relativisitic effects. A computation similar to 
that of Ebisawa et al. (1991), for Galactic black hole can- 
didates, has been done by Asaoka (1989) using the Kerr 
metric. However, our work incorporating both: the full gen- 
eral relativistic effects of rapid rotation, as well as realistic 



2 Bhattacharyya et al. 



equations of state describing neutron star interiors, using an 
appropriate metric is the first calculation for rotating neu- 
tron stars. In the present work we ignore the effects of the 
stellar magnetic field, so our results are applicable to weakly 
magnetised neutron stars. 

The structure of the paper is as follows. In section 2, 
we describe a framework for the calculation of neutron star 
structure, disc temperature profile and observed spectrum. 
We also comment here on the the chosen equations of state. 
We describe the numerical procedure for spectrum calcula- 
tion in section 3. The results of the spectrum calculation 
are displayed in section 4. We summarise the conclusions in 
section 5 and enlist some relevant mathematical expressions 
in the Appendix. 



2 THEORY 

2.1 Neutron steir structure calculation 

We calculate the structure of a rapidly rotating neutron star 

for realistic EOS models, using the same procedure as Cook 
et al. (1994) (see also Datta ct al 1998). With the assumption 
that the star is rigidly rotating and a perfect fluid, we take a 
stationary, axisymmetric, asymptotically flat and reflection- 
symmetric (about the equatorial plane) metric, given by 

dS"^ = g^.^dx'^dx'iiJL, = 0, 1, 2, 3) 

ojdtY 



-he ' '^r sm 



(1) 



where the metric potentials 7, p, a and the angular speed {uj) 
of zero-angular-momentum-observer (ZAMO) with respect 
to infinity, are all functions of the quasi-isotropic radial co- 
ordinate (f) and polar angle {0). The quantity f is related to 
the Schwarzschild-like radial coordinate (r) by the equation 
r = fe*-'^"''^^^. We use the geometric units c = G = 1 in this 
paper. 

We solve Einstein's field equations and the equation 
of hydrostatic equilibrium self-consistently and numerically 
from the centre of the star upto infinity to obtain 7, p, a, uj 
and f2 (angular speed of neutron star with respect to infin- 
ity) as functions of f and 8. This is done for a particular 
EOS and assumed values of central density and ratio of po- 
lar to equatorial radii. The obtained numerical equilibrium 
solutions of the metric enable us to calculate bulk structure 
parameters, such as gravitational mass (M), equatorial ra- 
dius (R), angular momentum (J) etc. of the neutron star. 
We also calculate the radius (rorb) of the innermost stable 
circular orbit (ISCO), specific energy (E) and specific angu- 
lar momentum (l) of a test particle in a Keplerian orbit and 
the Keplerian angular speed (Qk) (sec Thampan & Datta 
1998 for a detailed description of the method of calculation). 

2.2 Equation of state 

The neutron star structure parameters are quite sensitive to 
the chosen EOS. In the literature, there exist a large number 
of EOS ranging from very soft to very stiff. For the purpose 
of a general study, we have chosen four EOS, of which one 
is soft (EOS : A, Pandharipande (1979) (hyperons)), one 
is intermediate (EOS : B, Baldo, Bombaci & Burgio (1997) 
(AV14-h3bf)) in stiffness and two are stiff (EOS : C, Walecka 



(1974) and EOS : D, Sahu, Basu & Datta (1993)). Of these 
EOS : D is the stiffest. 

The structure of a neutron star for a given EOS is de- 
scribed uniquely by two parameters : the gravitational mass 
(A/) and the angular speed (Q). For each chosen EOS, we 
construct constant M equilibrium sequences with O varying 
from the non-rotating case upto the centrifugal mass shed 
limit (rotation rate at which inwardly directed gravitational 
forces are balanced by outwardly directed centrifugal forces). 
Depending on fi, M and the EOS model, neutron stars may 
have R > rorb or R < rorb- 



2.3 Disc temperature profile calculation 

The effective temperature of a thin blackbody disc is given 
by 



(2) 



where a is the Stefan-Boltzmann constant and F is the X- 
ray energy flux per unit surface area. We calculate F using 
the expression of Page & Thorne (1974) (valid for geomet- 
rically thin non-self-gravitating disc embedded in a general 
axisymmetric space-time) for the accretion disc around a 
rotating blackhole : 



where 



fir) 



(E - nKl)l,rdr 



(3) 



(4) 



Here nn is the disc inner edge radius and a comma followed 
by a variable as subscript to a quantity, represents a deriva- 
tive of the quantity with respect to the variable. 

Bhattacharyya et al. (2000) argue that, though formu- 
lated for the case of black holes, these expressions hold quite 
well also for the case of neutron stars. If 7? > rorb, then the 
disc will touch the star and we take rin — R. Otherwise, 

rin — rorb • 

The temperature profile, of the accretion disc calculated 
as described in this and the previous subsections is a func- 
tion of M and f2 of the central star for any adopted EOS. 



2.4 Calculation of the spectrum 

For calculating the disc spectrum in full general relativity 
including light bending effect for rapid rotation we adopt 
the following procedure. 

The observed spectrum from the accretion disc is given 

by 



F(Sob) 



(1/Sob) / 7ob(Kb)dnob 



(5) 



where the subscript 'ob' denotes the quantity in observer's 
frame, F is expressed in photons/sec/cm^/keV, E is photon 
energy in keV, I is specific intensity and 11 is the solid angle 
subtended by the source at the observer. 

As I/E'^ remains unchanged along the path of a photon 
(Misner et al. 1973), it is possible to calculate loh, if -^em is 
known (we use the subscript 'em' to denote the quantities 
in the emitter's frame). We assume the accretion disc to 
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radiate like a diluted blackbody (see Shimura & Takahara, 
1995). So Jem is given by 

/cm = {l/f)B{E,^,T,) (6) 

where / is the colour factor of the disc, B is Planck function 
and Tc = /Teff. 

The quantities -Eob and Eon are related by Ecm = 
/?ob(l + z), where (1 + z) contains the effects of both gravi- 
tational redshift and Doppler shift. (1 + z) is given by (Lu- 
minet 1979) 



1 + z = {1 + riKb sin a sin i){-gtt - '2QKgt,p 



-1/2 



(7) 



where i is the inclination angle of the disc with respect to 
an observer at infinity, b is the impact parameter of the pho- 
ton relative to the line joining the source and the observer 
and a is the polar angle of the position of the photon on 
the observer's photographic plate. Therefore bdbda is the 
apparent area of a disc element at observer's sky and the 
corresponding solid angle is given by 

b db da 

dn^b = — ^ (8) 

where D is the distance of the source from the observer. 
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Figure 1. Effect of general relativity: spectra from accretion disc 
around a neutron star of mass 1.4Mq. All the curves are for EOS 



model (B), Q = 0, D = 5 kpc. 



60°, M = 10^** g/sec and 



f = 2. Curve (1) corresponds to the Newtonian case, curve (2) to 
the general relativistic case including the effect of light bending 
and curve (3) to the general relativistic case without considering 
the effect of light bending. 



3 NUMERICAL PROCEDURE 

For a configuration, described by M and (and thus spec- 
ified by a set of Qiiu), we obtain Qk- To calculate the spec- 
trum for a given value of i with light bending effects, we 
backtrack the photons path from the observer to the disc, 
using standard ray tracing techniques (e.g. Chandrasekhar 
1983) and the relevant boundary conditions. For the met- 
ric nn) that we use, the equations of motion for photons 
are provided in the Appendix. We cover the disc between 
radii rin and rmid ~ lOOOrg; rin being the radius of the 
inner edge of the disc and rg the Schwarzschild radius (in- 
creasing Tmid has no significant effect on the spectrum) . Be- 
yond rmid, we ignore the effect of light-bending i.e., we take 
b sin a = r sin 4> ((/> is the azimuthal angle on disc plane) and 
dllob = {r dr d(j) cosi) / (Bhattacharyya et al. 2001). 

We have performed several consistency checks on our 
results: (1) by switching off the light bending effect (i.e. by 
considering flat space-time while backtracking the photon's 
path), we see that the spectrum matches very well with that 
computed by ignoring light bending effects (calculated by an 
independent code - Bhattacharyya et al. 2001). Also, in this 
case, the analytically calculated values of several quantities 
on the disc plane (e.g. r, (p, dcp/dt, dO/dt etc.) are reproduced 
satisfactorily by our numerical method, (2) an increase in 
the number of grid points on the (6,q) plane do not have 
any significant effect on the computed spectrum, (3) the 
spectrum matches very well with the Newtonian spectrum 
(Mitsuda et al. 1984) at low energy limit. This would imply 
that for higher frequencies, our spectrum is correct to within 
0.2% to 0.3%. 



4 RESULTS 

We calculate the general relativistic spectrum from the ac- 
cretion disc around rapidly rotating neutron star, taking into 



account the light bending effect. The spectrum is calculated 
as a function of 6 parameters : M, fi, distance of the source 
{D), inclination angle (i) (for face-on, i — 0°), accretion 
rate (Af) and colour factor /, for each of the chosen EOS. 
Our results are displayed in Figs. 1 to 5. In all the displayed 
spectra, we have assumed M = 1.4Mq (canonical mass for 
neutron stars), D = 5 kpc and / = 2.0. 

In Fig. 1, we have plotted the Newtonian spectrum and 
GR spectra with (LBGR) and without (NLBGR) light bend- 
ing effect, keeping the values of all the parameters same. At 
10 keV, the Newtonian flux is almost 2.5 times the LBGR 
flux. This is quite expected, because in the inner parts of the 
disc, Newtonian temperature is considerably higher than the 
GR temperature (see Fig. 2 of Bhattacharyya et al., 2000). 
LBGR flux is about 50% higher than NLBGR flux at 10 keV. 
This is because light bending causes the disc to subtend a 
larger solid angle at the observer than otherwise. Thus the 
general effect of light bending is to increase the observed 
ffux. 

According to Shimura & Takahara (1995), the thin 
blackbody description of the accretion disc, as adopted 
in this paper, is valid for O.IMe < M < MEdd, where 
Me = J/Edd/c^. Here I/Edd is the Eddington luminosity and 
A^Edd is the true Eddington accretion rate. For the purpose 
of demonstration, we have taken three different values of M 
in this range (for the mass shed conflguration) and plotted 
the corresponding spectra in Fig. 2. As is expected, we see 
that the high energy part of the spectrum is more sensitive 
to the value of M. It is seen that the spectra for different 
values of M are easily distinguishable. 

The inclination angle i is a very important parameter in 
determining the shape of the spectrum and its overall nor- 
malisation. In Fig. 3, we have plotted the spectra for four 
inclination angles, for the mass shed conflguration. We see 
that the observed flux at low energies is higher for lower 
values of i. This is due to the projection effect (proportional 
to cosi). But at higher energies (> 10 keV) this trend is 
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Figure 2. Accretion rate dependence; general relativistic spectra 
including light bending effects from accretion disc around a neu- 
tron star of mass shed limit configuration (Q = 7001 rad/s). Curve 
(1) corresponds to M = 10^^ g/sec, curve (2) to M = IQi^ g/sec 
and curve (3) to M = 2 X 10^^ g/sec. The values of all the other 
parameters are as in Fig. 1. 
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Figure 3. Inclination angle dependence: general relativistic spec- 
tra including light bending effects from accretion disc around a 
neutron star of mass shed limit configuration (Q = 7001 rad/s). 
Curve (1) corresponds to j = 0°, curve (2) to i = 30°, curve (3) 
to j = 60° and curve (4) to i = 85°. The values of all the other 
parameters are as in Fig. 1. 

reversed mainly because Doppler effect becomes important. 
The most energetic photons largely come from the inner 
portion of the disc, where the linear speed of accreted mat- 
ter is comparable to the speed of light. The net effect of 
Doppler broadening is a net blue shift of the spectrum, as 
a larger amount of flux comes from the blue-shifted regions 
than from the red-shifted regions. This is a monotonic trend, 
but it will be noticed from Fig. 3 that the curve for i = 85° 
overcomes that for i — 60° only at the edge of the figure, i. e. 
at energies > 30 keV. This is due to the fact that between 
these two inclinations the difference in the cos i factor is se- 
vere, and the blueshift overcomes this only at high energies. 

In Fig. 4, we have four panels for four inclination angles. 
In each panel, we have shown spectra for 3 different il (cor- 
responding to non-rotating, intermediate and the mass shed 
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Figure 4. Rotation rate dependence: general relativistic spectra 
including light bending effects from accretion disc around a neu- 
tron star. Panel (a) corresponds to i = 0°, panel (b) to i = 30°, 
panel (c) to i = 60° and panel (d) to i = 85°. In each panel, 
the solid curve corresponds to f2 = rad/s, the adjacent dashed 
curve corresponds to n = 7001 rad/s (the mass shed limit) and 
the dotted curve corresponds to H = 3647 rad/s. The values of 
all the other parameters are as in Fig. 1. 



configurations). With the increase of Q, disc temperature 
profile does not vary monotonically (see Fig. 3a of Bhat- 
tacharyya et al., 2000). Hence the behaviour of the spec- 
trum is also non-monotonic with SI. For non-rotating and 
mass shed configurations (for the assumed values of other 
parameters) the temperature profiles are very similar. As a 
result, the plotted spectra for these two cases lie almost on 
top of each other. However, for i = 0° the fiux corresponding 
to the mass shed configuration is slightly higher than that 
for f2 = 0, while the case is opposite at higher inclinations. 
This is a result of the inclination dependence of the (1 -I- z) 
factor given in Equation (^). 

In Fig. 5, we have compared the spectra for the four 
EOS models adopted by us, for configurations at the respec- 
tive mass-shed limits (which correspond to different values 
of Q because of the EOS dependence of the stellar struc- 
ture). The values of all other parameters have been kept 
the same. We see that the total flux received varies mono- 
tonically with the stiffness parameter, and is higher for the 
softer EOS. This effect has been noticed earlier by Bhat- 
tacharyya et al. (2000). We see that at high energies the 
fluxes for different EOS are considerably different. There- 
fore, fitting the observed spectra of LMXBs with our model 
spectra, particularly in hard X-rays, may provide a way to 
constrain neutron star EOS. Of course it must be remem- 
bered that these computations have been made assuming 
that the magnetic fleld of the compact object does not limit 
the inner boundary of the accretion disc. In the presence of a 
sufficiently strong magnetic fields, appropriate modifications 
must be taken into account while calculating the expected 
flux at high energies. 
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Figure 5. EOS model dependence: general relativistic spectra 
including light bending effects from accretion disc around a neu- 
tron star of mass shed configuration. Curve (1) corresponds to 
the EOS model (A) (H = 11026 rad/s), curve (2) to the EOS 
model (B) (H = 7001 rad/s), curve (3) to the EOS model (C) 
(Q = 6085 rad/s) and curve (4) to the EOS model (D) (Q = 4652 
rad/s). The values of all the other parameters are as in Fig. 1. 



5 CONCLUSION 

In this paper we have computed the observed radiation spec- 
trum from accretion discs around rapidly rotating neutron 
stars using fully general relativistic disc models. This is the 
first time such a calculation has been made in an exact way, 
without making any approximation in the treatment of ei- 
ther rotation or general relativity. In computing the observed 
spectrum from the disc, we explicitly include the effects of 
Doppler shift, gravitational redshift and light-bending for 
an appropriate metric describing space-time around rapidly 
rotating neutron stars. We find that the effect of light bend- 
ing is most important in the high-energy (> 3 keV) part 
of the observed spectrum. Photons at these high energies 
originate close to the central star, and hence their trajecto- 
ries are most affected by the light-bending effect. Depending 
on the viewing angle, this can enhance the observed flux at 
~ 10 keV by as much as 250% compared to that expected if 
light-bending effects are neglected. 

The calculations presented here deal only with the mul- 
ticolour blackbody disc. In reality, of course, there will be 
additional contributions to the observed spectrum from the 
boundary layer and also possibly from an accretion disc 
corona, both of which arc likely to add a power-law compo- 
nent at high energies (Popham & Sunyaov 2000, Dove et al 
1997). On the other hand, the spectra presented in Figs (2), 
(3) and (5) should remain essentially unaffected by bound- 
ary layer contribution, eis these are for neutron stars rotat- 
ing near the mass shed limit for which the boundary layer 
luminosity will bo negligible. For slowly rotating neutron 
stars, the disc component of the spectrum can be obtained 
by fitting and removing the contribution of the boundary 
layer, provided a good model for the boundary layer spec- 
trum is available. Popham & Sunyacv (2000) have calculated 
the boundary layer spectrum in the Newtonian approxima- 
tion. General relativistic modifications need to be included 
in these calculations to get a realistic estimate of the spec- 



trum of the boundary layer. We plan to address this issue 
in a future work. In the slow rotation case, the spectrum of 
the disc itself may be somewhat modified by the presence of 
a boundary layer if it extends beyond the disc inner radius 
assumed in our computations here, thus curtailing the inner 
edge of the disc. 

In addition to the contribution of the boundary layer, 
the possible contribution of an accretion disc corona to the 
emergent spectrum could also be significant. To be able to 
constrain the EOS models of neutron stars using the ob- 
served spectrum, this contribution must also be accurately 
estimated. In the present work, we restrict ourselves to thin 
black body and non-magnetic accretion discs in order to 
understand the effect of the EOS models describing neutron 
stars on the spectrum of the accretion disc alone (and thus 
neglect the effect of a corona) . We view this as the first step 
in modelling of spectra of accreting neutron stars including 
the effects of general relativity and rotation. 

The comparison of the non-rotating limit of our results 
with those of the fitting routine GRAD in the X-ray spec- 
tral reduction package XSPEC (Ebisawa et al. 1991), shows 
that the latter model overpredicts the high-energy compo- 
nent of the flux by a large factor. With the help of Ebisawa 
& Hanawa {private communication) we have been able to 
trace this disagreement to certain simplifying approxima- 
tions made in the GRAD code, as well as a couple of incor- 
rect expressions being used there. Conclusions based on the 
use of the GRAD routine may therefore need to be revised 
in the light of the new calculations presented here. 

The computation of the complete spectrum in the man- 
ner presented here is rather time-consuming and therefore 
not quite suited to routine use. Therefore, in order to make 
our results available for routine spectral fitting work, we in- 
tend to present a series of approximate parametric fits to 
our computed spectra in a forthcoming publication. 

The spectra presented here will find use in constraining 
the combined parameter set of the mass, the rotation speed 
and, possibly, the EOS, particularly of weakly magnetised, 
rapidly rotating neutron staxs. The relevant signatures are 
most prominent in hard X-rays, above ~ 10 keV. Sensitive 
observations of hard- X-ray spectra of LMXBs, therefore, are 
needed to fully utilise the potential of these results. 
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+2uJ{^l/{l - ^^))){d(j)/d\)(dt/d\) 
+e"-"''"(l-M')'((l/2)(7,.-p,.) 

~^i/{l~^^)){dcp/d\f 
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where, A is the affine parameter, L is the negative of the 
ratio of the (/i-component of the angular momentum and 
the f-component of the angular momentum of photon and 
a comma followed by a variable as subscript to a quantity, 
represents a derivative of the quantity with respect to the 
variable. 



APPENDIX 

For convenience, we use (= cos 9) instead of 9 and s (= 
f/{A + f)) instead of r as the coordinates. Consequently, the 
equation (|l|), becomes 

dS^ = -e^+''dt' + c^"((ylV(l-s)^)ds2 
+A^{s/{l^s)f{l/{l-^,^))df,^) 
+e^-M2(s/(l - s)f{l - ii'){d(l) - ujdtf (9) 

Here, A is a known constant of the dimension of distance. 

Now it is quite straight forward to calculate the geodesic 
equations for photons, which are given below. 

dt/d\ = e-<^+'''(l-(.^L) (10) 

+L/{e-'-^A\.s/{l-.s)f{l-f,')) (11) 

{ds/dXf = e-'"((l-s)VA')(e~'^+''>(l-c^L)^ 

-LV(c""''^'(s/(1-s))'(1-m'))) 
-s'{l-sf{l/{l-^'))y' (12) 

dfi/d\ = y (13) 

dy/dX = ~2{a,s + {l/{s{l~s))))y{ds/dX) 

+a,^(l/(s(l - s))f{l ~ ti^){ds/d\f 
-(a,M + (At/(1 - M^))y^ 



